Simulation of quantum dead-layers in ferroelectric tunnel junctions 
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In this letter, we simulate electronic transport through a metal-ferroelectric-metal tunnel junction 
by use of a nonequilibrium Green's function approach. We show that quantum effects such as Friedel 
oscillations lead to deviations from the Thomas-Fermi screening model. As a consequence, we predict 
a bistable resistive switching effect, depending on the polarization state of the ferroelectric tunnel 
barrier. 



The finite penetration depth of the electric field in a 
metallic system leads to considerable corrections to the 
electronic properties of a metal-insulator-metal tunnel- 
ing structure [jj]. For example, Mead |2| observed a de- 
viation from the classical geometric capacitance arising 
from the electric field penetration into the metal elec- 
trodes. This finding was explained theoretically by Ku 
and Ullman 0], as well as Simmons 0, by modeling 
the total capacitance as a sum of a (geometric) dielec- 
tric capacitance and electrode capacitance contributions 
(due to the finite field penetration). Black and Welser 
extended this approach to high-k materials (for exam- 
ple BaSr^Tii^Oa) In general, the field penetration 
effect is pronounced in capacitors and tunnel junctions 
containing ultrathin insulator films with a high dielectric 
permittivity. In addition, field penetration was also no- 
ticed in the context of finite size effects of complex oxide 
ferroelectric capacitors 0, H El ■ The finite screening 
length furthermore leads to uncompensated ferroelectric 
charge layers at metal-ferroelectric interfaces under short 
circuit conditions. First-principles calculations by Jun- 
quera and Ghosez predicted ferroelectric instabilities in 
SrRuOa/BaTiOa/Si'RuOa heterostructures with a tunnel 
barrier thickness even down to 2.4nm. (In fact, recent ex- 
periments on Lead compound oxides demonstrated ferro- 
electricity in ultrathin films [lfllllj |). In examples, the 
field penetration in the metal electrodes was described on 
the basis of the Thomas- Fermi screening length Q] . 

In this letter, we show that quantum interference ef- 
fects such as Friedel oscillations [l| lead to deviations 
from the semiclassical Fermi- Thomas screening within 
a metal-insulator-metal tunnel junction. Recently, we 
presented the concept of a ferroelectric tunnel junc- 
tion (FTJ) m 

based on a metal-ferroelectric-metal 
(MFM) layer sequence, consisting of two (semi-) metals 
separated by an ultrathin ferroelectric layer as visualized 
in the inset of Fig.QJa). 

In order to simulate the IV characteristics of the FTJ, 
a no n-eq uilibrium Green's function approach has been 
used \lA I If* 1 1 (J ] . The employed empirical tight-binding 
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FIG. 1: Simulated electron density distribution (a) and po- 
tential profile (b) (with LDOS grayscale plot) for zero bias. 
(The Fermi energy corresponds to the dash-dotted line.) The 
inset shows a schematic sketch of the tunnel junction. 



one-band formulation constitutes a good approximation 
for the description of conduction band electrons inside 
the ferroelectric tunnel barrier and the (semi-) metal- 
lic electrodes. Within this framework, Coulomb charg- 
ing effects are accounted for in terms of a self-consistent 
Hartree-potential. Polarization charges at the two in- 
terfaces of the ferroelectric material are implemented as 
charged mono-layers. (For simplicity, we assume a mono- 
domain polarization state with a polarization vector per- 
pendicular to the junction plane. The device structure 
is translationally invariant within the grown layers.) An 
applied bias voltage is implemented as a difference in the 
local Fermi-energies of the two reservoirs. 

In the following, we will discuss an idealized case of 
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an MFM structure, assuming an effective electron mass 
of m* = 0.5m e and a dielectric constant of e r = 300 
throughout the device. The ferroelectric tunnel barrier 
of 3.2nm and 0.5eV height (with interface polarization 
charges of a = ±5 x 10 13 cm~ 2 ) is surrounded by two 
(semi-) metallic contact regions with electron concen- 
trations of 10 21 cm -3 and 10 20 cm~ 3 , respectively. (Of 
course, the chosen parameters have to be adjusted to the 
individual experimental device structure.) 

Fig.^a) shows the simulated electron contration under 
zero bias conditions (short circuit case). Due to the wave- 
nature of the electron, Friedel oscillations arise within 
the electrode regions in the vicinity of the tunnel barrier 
(clearly pronounced in the left contact). The oscillation 
period is given by the Fermi wavelength divided by two 
(Af/2 w 1.0 nm in the l.h.s. reservoir) |y. As an im- 
portant consequence, quantum mechanically determined 
partially depleted regions close to the tunnel barrier arise 
(see arrows in Fig. Ufa): l.Onm and 2.2nm), in contrast 
to the abrupt electron concentration drop in the classical 
case. We refer to this phenomenon as a "quantum dead- 
layer" (analogous to a conventional dead-layer which 
also introduces a non-screening and non-ferroelectric re- 
gion at the metal- insulator interface). Within the semi- 
metallic electrode on the r.h.s., the depleted region is 
further extended due to band-bending, varying with the 
applied voltage. In Fig. IHb), the corresponding self- 
consistent potential profile is shown (solid line) with a 
grayscale plot for the local density of states (LDOS). 
One can clearly identify the polarization charge at the 
interfaces of the ferroelectric tunnel barrier in terms of a 
kink in the continuous part of the potential and a tilted 
barrier potential due to the resulting depolarization field. 
The LDOS plot reveals the wave-properties of the elec- 
tron: Firstly, standing wave patterns result from reflec- 
tions at the tunnel barrier, giving rise to the Friedel os- 
cillations with a quantum dead-layer. Secondly, multiple 
reflections within the barrier region yield resonances in 
the higher energy region. (Latter may become relevant 
for hot-carrier injection conditions.) Furthermore, within 
the reservoir on the l.h.s., a shallow triangular quantum 
well arises, which is responsible the screening charge. 

The local electric field within the ferroelectric barrier 
as a function of the applied voltage is shown in Fig.Efa). 
Here, an intrinsic coercitive field of 700kV/cm has been 
assumed, resulting in transition voltages of approxi- 
mately +0.6V and —0.5V. The asymmetry stems from 
the assumed electron concentration asymmetry within 
the electrodes. Since there is a significant voltage drop 
within the electrode regions close to the barrier (depleted 
regions), the external voltage necessary to induce a po- 
larization flip in the ferroelectric material is significantly 
higher than the intrinsic voltage drop across the barrier. 
Additionally, in Fig. El[b), the corresponding current- 
density vs. voltage characteristics is plotted in a log- 
scale, revealinge an exponential behavior of the current 
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FIG. 2: Simulated electric field (a) and current-density 
(b) characteristics. (The assumed intrinsic coercitive field is 
700kV/cm.) 

density which is determined by the tunnel barrier. De- 
pending on the the polarization state and the sign of the 
applied voltage, the slope in the log-scale varies due to 
different potential profiles. The bistable IV characteris- 
tics could make such an MFM device a possible candidate 
for a non- volatile memory element. Obviously, the total 
capacity of the considered device structure will also be 
affected by the quantum dead-layer and depletion effects 
discussed above. 

In conclusion, we have considered quantum interfer- 
ence effects in ultrathin metal-ferroelectric-metal tunnel 
junctions. The current-voltage charateristics of an exem- 
plary tunnel structure has been simulated with the help 
of a noncquilibrium Green's function formalism. As a 
result, a bistable resistive switching behavior was pre- 
dicted. 
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